Rotating Fulde-Ferrell-Larkin-Ovchinnikov state in cold Fermi gases 
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We study an effect of rotation on the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state of two 
component Fermi superfluid gases in a toroidal trap. We investigate a stability of tlie FFLO states 
in the quasi-one-dimensional regime on the basis of the Bogoliubov-de Gennes equation. We find 
that two novel FFLO phases, i.e., the half quantum vortex state and the intermediate state of Fulde- 
Ferrell (FF) state and Larkin-Ovchinnikov (LO) state, are stabilized by the rotation. The phase 
diagram for the FF state, LO state, intermediate state, and half quantum vortex state is shown in 
both T-P plane and T-h plane. We demonstrate characteristic features of these states, such as the 
order parameter, fiux quantization, and local polarization. Several related works are discussed, and 
the advantages of cold Fermi gases are indicated. 

PACS numbers: 67.85.-d, 71.10.Ca, 74.25.Dw 
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I. INTRODUCTION 

Two component Fermi gases with population imbal- 
ance attract much attention from both theoretical and 
experimental point of view [H-Q. One of the moti- 
vations of current studies is the realization of Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state [1, Q, which 
is characterized by a spatial oscillation of the order pa- 
rameter arising from the center-of-mass momentum of 
Cooper pairs. This state is induced by a mismatch of 
the Fermi surfaces between two component Fermi par- 
ticles produced by a population imbalance. Since many 
parameters, such as the population imbalance, interac- 
tion strength, and trap potential, can be experimentally 
controlled 3] , cold Fermi gases are regarded as a promis- 
ing candidate for the FFLO superfluid. Recent experi- 
ment found an evidence for the FFLO superfluid state in 
the elongated harmonic trap , and attract much atten- 
tion in a variety of fields such as the condensed matter 
physics (rl-flU and nuclear physics [163 where possible 
FFLO phases have been discussed. 

The FFLO state has an internal degree of freedom aris- 
ing from the inversion symmetry, and therefore, there 
are several phases of the FFLO state. Most of the- 
oretical works on the FFLO state are focused on the 
two phases. One is the Fulde-Ferrell (FF) state which 
is described by the order parameter A(r) oc e^'^ '", and 
the other is the Larkin-Ovchinnikov (LO) state in which 
A(r) oc cos(q • r). The LO state is regarded as a mixture 
of two FF states with opposite momentum A(r) cx e*'' '' 
and A(r) oc e~"i "". It has been shown that the LO state 
is stable against the FF state in most cases [§]. 

The FFLO state in cold Fermi gases has been studied 
in the literature [l7l - [23 |. Some of their works investi- 
gated the effect of harmonic trap potential which plays an 
essential role for the spatial structure of the FFLO phase. 
They found that the radial FFLO (R-FFLO) state is sta- 
ble, in which the order parameter changes its sign along 
the radial direction around the edge of the harmonic trap. 
The spatial symmetry remaining in the harmonic trap, 
however, is not broken in the R-FFLO state, thus it may 



be difficult to experimentally observe the spatial oscilla- 
tion of order parameter. On the other hand, one of us 
showed that the angular FFLO (A-FFLO) state in which 
the order parameter changes its sign along the angular di- 
rection is stable in a toroidal trap 25] . Since the A-FFLO 
state breaks the rotation symmetry, one expects that a 
spatial modulation characteristic of the FFLO state can 
be detected in experiments. Furthermore, A-FFLO state 
is a novel FFLO state in the sense that the rotation sym- 
metry is spontaneously broken. This is in sharp contrast 
to the superconductors where the rotation symmetry is 
broken by the crystal lattice and therefore the A-FFLO 
state is not stabilized. Hence, it is interesting to study 
the response of A-FFLO state to the rotation. It is ex- 
pected that another novel FFLO state is stabilized by the 
rotation. A related work has been given by Kashimura et 
al. They showed that a spontaneous current appears in a 
toroidal trap with a potential barrier [26| , but they didn't 
investigate the effect of rotation. In this paper, we inves- 
tigate the FFLO phases induced by the rotation. It is 
expected that those studies will be tested by the experi- 
ments since such controllable experiment is an advantage 
of cold atom gases. 

An effect of rotation on the FFLO state has been stud- 
ied in the context of cold Fermi gases [27l - [29| as well as 
in mesoscopic superconductors where the ma gne tic field 
plays an equivalent role to the rotation [2^ Isol l3l| . It has 
been shown that the nucleation of vortex in the FFLO 
superfluid gives rise to intriguing phenomena. In con- 
trast to those works, our study is focused on the quasi- 
one-dimensional superfluid formed on the toroidal trap, 
which is an ideal system to study the effect of gauge fleld 
on the FFLO state. The A-FFLO state is regarded as 
the LO state along the quasi-one-dimensional ring, while 
order parameter of the vortex state is the same as that 
of the FF state. In this work, we study how the A-FFLO 
state formed on the toroidal trap is changed by the nu- 
cleation of vortex. The same issue has been investigated 
in the context of mesoscopic superconducting ring un- 
der the tilted magnetic fleld 132!], but their study over- 
looked many important results, as we will show later. 
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We show that the "giant vortex state" corresponding to 
the FF state is stabiHzed in the imbalanced gas near the 
superfluid critical temperature Tc, and the phase transi- 
tion to the "intermediate state" between the FF and LO 
states occurs with decreasing the temperature. Further- 
more, the "half quantum vortex state" becomes stable in 
a certain range of the rotation and population imbalance. 
This phase leads to the half quantized flux of mass and 
the unusual Little-Parks oscillation of Tc. We elucidate 
superfluid properties of these FFLO states, such as the 
superfluid order parameter, local polarization, and the 
flux of mass. 



II. FORMULATION 

We study rotating two component Fermi gases in a 
toroidal trap, in which atoms are loaded on a quasi-one- 
dimensional ring [s^, [sj] ■ We assume that the angular 
velocity is perpendicular to the plane of ring, leading to 
O = fiz. Then, the Hamiltonian is approximately de- 
scribed by the following one-dimensional attractive Hub- 
bard Hamiltonian: 
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(1) 
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'JO- — Cj^Cja 



where a =t, i describe two hyperfine states, n 
is the number operator, U is the coupling constant of 
attractive interaction, and /Xo- is the chemical potential 
for cr particles. The mass of atoms M is related to t — 
1/2M, and the phase $ = Rfl/2t denotes the Peierls 
phase, where R is the radius of ring (see Appendix A). 
We take the unit ?i = fee = 1 and t = 1. 

Although one-dimensional models can be solved with 
use of sophisticated analytic or numerical methods (ssl - 
l42j . we adopt the Bogohubov-de Gennes (BdG) equa- 
tion, because our interests are focused on the quasi-one- 
dimensional regime where a three-dimensional long range 
order is achieved [isUisj . Then, the Hamiltonian is re- 
duced to the mean field Hamiltonian 

Hra = -I- e~'*c]^Cj + i^) - H^h^a 



u 



(2) 



where Aj = —\U\{cj]^Cj-^) is the order parameter. We fo- 
cus on the weak coupling BCS region, and then the BdG 
equation is quantitatively appropriate. The Hartree term 
can be neglected there. We don't touch the BCS-BEC 
crossover region where higher order corrections beyond 
the mean field theory play important roles d, |47| , 
although it is expected that the FFLO state is robust 
there l25ll48t. 



The mean field Hamiltonian is diagonalized with use 
of the following Bogoliubov transformation: 



(3) 



where 7^ and 7,^ are the creation and annihilation opera- 
tors of quasiparticles, respectively. The wave function of 
quasiparticles satisfies the BdG equation 
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where Hji„ = -te**5j_i,i - te **(5j+i,/ - ^icrSji- 

The self-consistent equation for the order parameter is 
obtained as. 



A, = -\U\ 



(5) 



where f{E^) is the Fermi distribution function. We can 
refer to the analytic solution of BdG equation in one di- 
mension [4^, [13], however, we numerically solve the self- 
consistent equation (Eqs. @ and ([S])) to investigate the 
thermodynamic stability of superfluid phases. Several 
self-consistent solutions corresponding to the metastable 
states are obtained. We determine the stable phase by 
comparing the free energy of those states. The free en- 
ergy is evaluated to be 



F = E-TS 



|A,| 



\U\ 



+TY{f(-E^)^^s[fi~E,)]+fiE,)log[fiE,)]}. 

(6) 

We solve BdG equations for a fixed chemical poten- 
tial pL„ in the ground canonical ensemble and calculate 
the particle number and population imbalance which are 
expressed as N = N-^ + and P = (iV^ - N^)/N, re- 
spectively. We obtain the number of a particles and 
iVj. as 



^t = E^^t= EE 
3 3 V 

^i = E-.i = EEi 
3 3 V 



"it I 



1 - !{EJ)\ 



(7) 



respectively. We define the "magnetic field" h = {fj,^ — 
/i4,)/2 as the difference of chemical potential between t 
and 4 particles. 



III. NUMERICAL RESULTS 

We set \U\/t = 1.5 and assume the chemical potential 
= (/x-f -I- = -0.8) leading to rijcr = (fija) < 0.4 

so that the lattice Hamiltonian appropriately describes 
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the continuum gas. The fohowing results are obtained 
for the number of lattice sites Nj^ = 200 with imposing 
the periodic boundary condition. Note that the periodic 
boundary condition should be satisfied for a gas on the 
ring. 

In all our results the order parameter Aj is approxi- 
mately described as 



A, 



A_e-*«-^ 



(8) 



where the center-of-mass momentum of Cooper pairs is 
q± = 2i:m±/Ni^ with m± being an integer so as to sat- 
isfy the periodic boundary condition. The FF state is 
described by a set of order parameter (A^., A_) oc (1,0) 
or (A+,A_) oc (0,1), while (A+,A_) oc (1,±1) with 
TO-i- = m_ = TO in the LO state. 



A. Rotating FFLO phases 

Before studying the effect of rotation, we discuss the 
phase diagram of gases at rest. In Fig. [1] we show the 
phase diagram for $ = in the T — h plane as well 
as in the T — P plane. We denote the superfluid crit- 
ical temperature Tc{h,^/^o) as a function of the mag- 
netic field h and the normalized phase factor "I>/<I>0: where 
$0 = 27r/A^L- The temperature in Fig. 1 is normalized 
by the critical temperature of the balanced gas at rest, 
Tc(0,0) = 0.0816. Figuredla) shows that the BCS state 
is stable for h < 0.0915 while the LO state is stabilized 
for h > 0.0915 at T = 0. Note that the order param- 
eter of the LO state is described as A^- a cos[{^om)j]. 
Since we consider the Fermi gas trapped on the ring, this 
LO state corresponds to the A-FFLO state studied in 
Ref. [1^. The integer to increases with increasing the 
magnetic field h [9]. Although the FFLO state is stable 
for a huge magnetic field h ^ Tc(0, 0) in our calculation, 
our study focuses on the moderate magnetic field h ^ Tc 
where to <C A^l- Figure 1(b) shows that the BCS state 
is unstable in the imbalanced gas P > at T = 0. This 
is because the BCS state does not have an excess parti- 
cle owing to the excitation gap. We also find forbidden 
regimes (FR) in the T — P plane since the particle num- 
ber is finite and not fixed in our calculation based on 
the grand canonical ensemble. Another inhomogeneous 
superfluid state may be stabilized in the canonical en- 
semble, but that is beyond the scope of our paper. 

We here turn to the rotating systems. Figure ^a.) 
shows the phase diagram in the T — $ plane for a field 
h = 0.13 where the LO state with to = 4 is stable at rest, 
while Fig. [21b) shows the phase diagram for h = 0. While 
a conventional vortex state and usual Little Parks oscil- 
lation of Tc appear in the rotating BCS state (Fig.[2l^b)), 
we see several intriguing features in the rotating FFLO 
state (Fig. [^a)). First, the FF state with A+ = or 
A_ = is stabilized near the critical temperature in the 
rotating gas with $/$o 7^ 0. The two FF states A^ oc e**?^ 
and Aj oc e"*''-' are degenerate at rest, and therefore, the 
mixture of them, i.e., the LO state Aj oc (e*^-' -t-e"*"^^ ) is 



stabilized at T = Tc. On the other hand, the degeneracy 
is lifted by the rotation, and then one of the FF states 
is stabilized near T^. This FF state is regarded as a "gi- 
ant vortex state" because the vorticity to is much larger 
than the conventional one. The critical temperature is 
increased by the rotation since one of the FF states is 
stabilized. The FF state may be furthermore stabilized 
by the Hartree term, which is neglected here, because of 
the Fermi liquid correction [sij . 

Second, the "half quantum vortex state" is stabilized 
by the rotation. Figure[2](a) shows that the LO state with 
(to+,to_) ~ (4,4) changes to the phase (to+,to_) — 
(5,4), as the rotation, i.e., the phase $ increases. The 
latter phase is regarded as the half quantum vortex state 
since the order parameter is approximately described at 
low temperatures as 



Aj oc exp[i($o X 5)j] exp[-i($o x 4)j] 
oc exp[i($o X 0.5) j] cos[(<I>o x 4.5)j]. 



(9) 



The exponential provides the phase tt for a circuit along 
the ring, which corresponds to a half quantized vortex. 
Another phase tt required for the uniqueness of order pa- 
rameter is given by the odd number of phase changes in 
the sinusoidal function. Then, the order parameter of the 
half quantum vortex state is regarded as the product of 
those in the FF and LO states with an unconventional 
period. Similarly, the phase with (to+,to„) = (6,3) 
is the half quantum vortex state while the phases with 
(to_(-,to_) — (5,3) and {m^,m^) = (6,2) are the LO 
state having one and two quantized vortex, respectively. 
Thus, the half quantum vortex state is realized when the 
difference of vorticity is odd between the coexisting FF 
states. One of the two quantum numbers to+ and to_ 
changes with the increase of rotation, and induces the 
half quantized vortex. This is in sharp contrast to the 
rotating BCS state which has only one quantum number 
TO. In other words, the half quantized vortex is induced 
by the multi-component order parameters of the FFLO 
state, as in the chiral p-wave superconductor [52| . 

We see a characteristic feature of the half quantum 
vortex state appearing in the quantized flux of mass. We 
show the rotation dependence of the flux of mass / for 
h ^ and = 0.13 at T = in Fig. [3l which is obtained 
as 



f = h + A, 



= 2tQ 



(10) 

(11) 



/, = 2t3|^u^,<;i,[l - f{E,)]\. (12) 

Since the flux is independent of the site, we omit the 
index j in Eqs. ([TT]) and (|12l) . Figure [3| clearly shows 
that the flux of mass is half quantized in the rotating 
FFLO state {h = 0.13). 

Owing to the nucleation of half quantum vortex, the 
shows an unconventional Little-Parks oscillation studied 
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FIG. 1. (a) T — h and (b) T — P phase diagram for the gas at rest (^/^o = 0). The integer m denotes a center-of-mass 
momentum of Cooper pairs q — <l>om. The solid and dashed lines show the second and first order transition lines, respectively. 
The BCS superfluid state with m = is denoted as "BCS" , while m > for the LO state is shown in the figure. The character 
"FR" denotes the forbidden regimes in Fig. 1(b). The LO state with m = 1 is stabilized in a tiny region of the T — h plane, 
but it cannot be resolved in Fig. 1(a). 
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FIG. 2. r — $ phase diagram at (a) h = 0.13 and (b) h = 0. The solid and dashed lines show the second and first order 
transition lines, respectively. The integers m± denote the momentum of Cooper pairs, and A± are the amplitude of order 
parameter for m±, respectively (see Eq. ([5])). The phase diagram is periodic for <I> with the period 3>o (see Appendix B). 



in Ref. [sTj . Figure [2][a) shows the Little Parks oscilla- 
tion with the period being a half of the conventional one. 
It should be noticed that the amplitude of oscillation is 
remarkably enhanced in the FFLO state compared with 
the BCS state. 

Third, we point out another phase in Fig. [2][a), that 
is the mixture of LO state and FF state. For instance, 
the FFLO state with (m+, m_) = (4, 4) is not a pure LO 
state in the sense that the amplitude A_|_ for m+ — 4 
is not equivalent to A_ for to_ =4. The temperature 
dependence of the ratio S = A+/A_ at h = 0.13 and 
<i>/<i>o = 0.0637 is shown in Fig. SI where the superfluid 
phase is divided into three regions, (I), (II), and (III). 
In the low temperature region (I), the ratio 6 is almost 
unity corresponding to the LO state, while the FF state 
with (5 = is stabilized in the high temperature region 
(III). Interestingly, the intermediate ratio < (5 < 1 is 
obtained in a wide intermediate temperature region (II). 
Thus, the "intermediate state" between the FF and LO 
states is stabilized by the rotation, and then the order 



parameter is described as 

Aj oc (5exp[i($o X 4)j] + exp[-i($o x 4)j]. (13) 

Although the intermediate state changes to the LO 
state through the crossover, these states can be distin- 
guished by observing the flux of mass. As shown in 
Fig. 21 the flux of mass is almost zero in the LO state. 
This is because the susceptibility to the rotation is sup- 
pressed by the excitation gap. A narrow band is formed 
by the Andreev bound states localized around the nodes 
of order parameter (45j . its contribution is negligible 
for these parameters. On the other hand, the flux of 
mass increases with temperature in the region (II), and 
shows a peak at the transition temperature to the FF 
state, T = Tc2 ^ 0.946Tc(0. 13, 0.0637). The flux de- 
creases as approaching to the superfluid critical temper- 
ature T = Tc because the order parameter of FF state 
is decreased. In other words, the flux of mass is first 
enhanced below owing to the formation of the giant 
vortex state, and then suppressed below Tc2 as decreas- 
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FIG. 3. (Color online) The rotation dependence of the flux of 
mass at r = in the BCS state {h = 0, blue solid line) and 
in the FFLO state {h = 0.13, red dotted line). The flux of 
mass is normalized by the value at /i = and $/$o = 0.5. 
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FIG. 4. Temperature dependence of the ratio 5 = A+/A-_ 
(solid line) and the flux of mass / (dashed line) at h — 0.13 
and ^/^o = 0.0637. 

ing the temperature. Such nonmonotonic temperature 
dependence is a characteristic feature of the FFLO state, 
while that is not observed in the BCS superfluid. Note 
that we here assume the system without conserving the 
angular momentum. Another phase diagram will be ob- 
tained when the angular momentum is conserved across 
the superfluid transition. 

B. Local polarization 

The FFLO phases studied in the previous subsection 
show their characteristic feature in the local polariza- 
tion. Figure [S] shows the spatial profile of local polariza- 
tion pj = rij-f — rij^ in the FF, LO, half quantum vor- 
tex, and intermediate states. It is shown that the local 
polarization in the LO state shows peaks at the quasi- 
nodal points of order parameter where | | ~ (compare 
Fig. E] with Fig. [6]), as in the A-FFLO state [23]. This is 
because excess particles are localized around the quasi- 
nodal points of order parameter in order to make the 
loss of condensation energy as small as possible. There- 



fore the number of peaks in the local polarization uq can 
be regarded as the number of quasi-nodal points in the 
order parameter. On the other hand, the FF state shows 
a uniform polarization, since the translation symmetry is 
not broken there. The spatial inhomogeneity is smoothly 
enhanced by decreasing the temperature in the interme- 
diate state. Thus, we can discriminate the LO state, in- 
termediate state, and FF state by the polarization mea- 
surements. 

A characteristic feature of the half quantum vortex 
state appears in the number of peaks in the local po- 
larization. The number uq is odd (even) in the half (in- 
teger) quantum vortex state. (See above discussion and 
Fig. [6]) For example. Fig. [5ja) shows uq = 8, while that 
is riQ = 9 in Fig. \^h) where the amplitude is described 
by \Aj\ cx I cos[($o x 4.5)j]|. Thus, the odd number of 
no is a clear signature of the half quantum vortex state. 
Note that the no is the same as the number of excess 
particles. As increasing the excess particles, the number 
of quasi-nodes increases so as to hold them. 



C. Mesoscopic effect 

We here comment on the magnetic field dependence 
of superfluid phases. Although the LO state with 
(m+,m_) = (4,4) is stabilized at rest in both magnetic 
fields h — 0.115 and h — 0.13, the FF state induced by 
the rotation depends on the field. While the FF state 
with A+ = is stabilized for h = 0.13 (Fig. [IJa)), the 
other FF state with A_ = is stabilized for h = 0.115 
as shown in Fig. [T] Accompanied by the change of FF 
state, the half quantum vortex state also changes as in- 
creasing the magnetic field from h = 0.115 to h = 0.13; 
a small rotation stabilizes the half quantum vortex state 
(m+, m_) = (4, 3) in the former, while (m+, m_) = (5, 4) 
in the latter. 

We understand these magnetic field dependences in the 
following two ways. First, the superfluid phase having a 
large number of zeros of order parameter no is stabilized 
by increasing the magnetic field so as to gain the mag- 
netic energy. For example, no = 7 [(m+,r7i_) — (4,3)] 
a.t h = 0.115 increases to no = 9 [(m+,m_) = (5,4)] at 
h = 0.13. 

Another interpretation is based on the discrete single 
particle spectrum, that is the mesoscopic effect. Figure[H] 
shows the schematic sketch of the Cooper pairing for the 
cases h = 0.115 and h — 0.13. The upper panel shows 
the Fermi momentum at rest, while the lower panel shows 
the shift of Fermi momentum due to the rotation (see Ap- 
pendix A) . This shift gives rise to the lifting of degeneracy 
of two FF states with A+ = and A_ = 0. In Fig.IHJa), 
a small shift moves the Cooper pairs with m_ — 4 away 
from the Fermi surface. Since the Cooper pairs should 
be formed by the quasiparticles near the Fermi surface, 
the FF state with m+ = 4 is favored at low magnetic 
fields as shown in Fig.[51[a), while the other FF state with 
m_ = 4 is stabilized at high magnetic fields as shown in 
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FIG. 5. (Color online) Spatial profile of the local polarization at h — 0.13. (a) ^/^o = 0.0637. The dotted (red), solid 
(green), and dashed (blue) lines represent the data for T/rc(0.13, 0) = (LO state), r/rc(0.13, 0) — 0.947 (intermediate state), 
and r/rc(0.13, 0) = 0.991 (FF state), respectively, (b) $/$o = 0.159. The dotted (red), solid (green), and dashed (blue) 
lines represent the data for r/rc(0.13,0) = (half quantum vortex state), r/rc(0.13, 0) — 0.969 (intermediate state), and 
T/Tc(0.13,0) = 1.01 (FF state), respectively. 
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FIG. 6. (Color online) Spatial profile of the amplitude of order parameter \Aj \ at h = 0.13. In (a) and (b), the parameters are 
the same as in Figs. [5]^ a) and (b), respectively. 
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FIG. 7. T — $ phase diagram at h — 0.115. The other pa- 
rameters are the same as Fig. [5] 



Fig.HJb). In other words, the reconstruction of quasipar- 
ticle energy due to the rotation determines the stable FF 
state near Tc. When the temperature is decreased, the 
sub-dominant order parameter appears below Tc2- That 
is A_|_ with rn+ = 4 (A_ with rn_ = 4) for a small rota- 



tion, and changes to m+ = 5 (rri- = 3) as increasing the 
rotation. Thus, the FF state with A+ = changes to 
the half quantized vortex state (m+,TO_) — (5,4), while 
the other FF state with A_ — changes to another half 
quantum vortex state (m+,m_) — (4,3). Note again 
that these differences of the T— $ phase diagram between 
h — 0.115 and h = 0.13 are induced by the mesoscopic 
effect. 



D. Phase diagram 

Finally, we discuss the phase diagram of the rotating 
Fermi superfluid gases in the T — h and T — P plane. We 
fix the rotation $/4>o = 0.286 in Fig. [9] For this param- 
eter the BCS state is changed to the vortex state with 
vorticity m — 1. More important finding is that the half 
quantum vortex state [for example, {m^,m^) = (1,0)] 
is induced in a certain range of the population imbal- 
ance. Therefore, the flux of mass iteratively changes as 
1 — i — > 1 — > i with increasing the population imbal- 
ance. Another point is that the FF state is stabilized near 
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FIG. 8. Schematic sketch of Cooper pairing in the FFLO state, (a) the case of ft = 0.115 and (b) the case of ft = 0.13. The 
upper panel shows the quasiparticle energy of the gases at rest, while the lower panel shows the shift of Fermi momentum 
due to the rotation. The circles represent the momentum discretized by the mesoscopic effect, k^a and k^^ are the Fermi 
momentum for a particles. 
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FIG. 9. (a) r — ft and (b) T — P phase diagram for $/<l>o = 0.286. Scaled-up figures of (a) and (b) are shown in (c) and (d), 
respectively. 



the critical temperature, as studied in section IIII Al In 
order to illuminate these features, we show the phase di- 
agram scaled-up in the large population imbalance (mag- 
netic field) region. We see that the population imbalance 
changes the vorticity on a sequence of phase transition 



lines. These features are understood on the basis of the 
mesoscopic effect discussed in section UlI CI 

For a rotation smaller than <i>/<f>o = 0.286, where the 
vortex is not nucleated in the balanced gas, a half quan- 
tum vortex can be induced by the population imbalance 
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(compare Fig.[2lja) with Fig. [2Kb)). Then, the gas begins 
to rotate by increasing the population imbalance, if the 
angular momentum is not conserved for some reasons. 
This intriguing phenomenon would be a signature of the 
FFLO superfiuid gas. 



IV. SUMMARY AND DISCUSSION 

We have studied effects of rotation on the FFLO super- 
fluid state of imbalanced Fermi gases loaded on a quasi- 
one-dimensional toroidal trap. We found several novel 
phase transitions induced by the rotation and population 
imbalance. First, the FF state is stabilized in the rotating 
gases near the critical temperature, although this state is 
unstable against the LO state in the gases at rest. This 
means that the "giant vortex" appears in the imbalanced 
gases near because the FF state on the ring is regarded 
as a vortex state. We found that the vorticity shows a 
nonmonotonic and intriguing dependence on the popu- 
lation imbalance. Second, the FF state changes to the 
intermediate state through the second order phase tran- 
sition at T = Tc2 < Tc. As decreasing the temperature 
below Tc2, the crossover occurs from the intermediate 
state to the LO state. We point out that the FF, LO, and 
intermediate phases are distinguished by the flux of mass 
and local polarization. Third, we showed a sequence of 
quantum phase transitions with increasing the rotation. 
The FFLO state with an integer vortex changes to the 
half quantum vortex state. Then, the flux of mass is half 
quantized (Fig.|3]). The half quantum vortex state arises 
from the multi-component order parameters of the FFLO 
state, and therefore the half quantized flux would be a 
clear evidence for the FFLO superfiuid state. The spatial 
profile of the order parameter in the half quantum vortex 
state is described by the product App x Alo where App 
and Alo are the order parameters of FF and LO states, 
respectively. 

These results are obtained on the basis of the BdG 
equation. Since the BdG equation is a mean field approx- 
imation it breaks down in the purely one-dimensional 
system. It is well known that the long-range order of su- 
perfluidity/superconductivity is fragile against the fluc- 
tuation in one dimension. We adopt the one-dimensional 
model for simplicity, however our interests are focused 
on the quasi-one-dimensional superfiuid which is realized 
in the toroidal trap [lH, [s^ 113] • It has been shown that 
a weak three-dimensional coupling gives rise to the long 
range order fist-lisj. The BdG equation is qualitatively 
appropriate for a weak attractive interaction in this case. 
The correlation effect beyond the mean field approxima- 
tion plays important roles for a strong attractive interac- 
tion. Therefore, our study basically focuses on the weak 
coupling BCS region. However, it has been shown that 
the FFLO superfiuid state is robust against the strong at- 
tractive interaction in the unitary regime [25l IZsj . This 
indicates that the novel FFLO phases studied in this pa- 
per are robust in the BCS-BEC crossover region. 



We note that a fingerprint of the FFLO phases appears 
in the purely one-dimensional system although the 
long range order is suppressed there. For instance, the 
free energy shows a characteristic periodicity with respect 
to the phase $ as shown in Ref . [5j| . We expect that the 
free energy of the imbalanced gases shows an oscillation 
with respect to the phase $, and the period is a half of 
the BCS state. This oscillation is a fingerprint of the 
half quantum vortex state in the FFLO state. It will be 
interesting to see this periodicity with use of the accurate 
numerical and/or analytical methods for one-dimensional 
models. 

Three comments are given on the related theoretical 
works. First, the same model as in Eq. ([T]) was solved in 
the recent paper [s^ on the basis of the BdG equations. 
However, they missed most of the novel FFLO phases 
induced by the rotation, except for the vortex state in 
the small imbalance region. We showed T — h and T — P 
phase diagrams in which the FF state with giant vortex, 
half quantum vortex state, and intermediate state are 
stabilized. These phases were not shown in Ref. [s^ . 

Second, a superfiuid state similar to the A-FFLO state 
was investigated in the optical lattice by Chen et al. |53] . 
The response to the rotation is quite different from our 
case, because the rotation symmetry is broken by the op- 
tical lattice potential. For example, the FF state and the 
intermediate state are suppressed because the degener- 
acy of the FF states is lifted. In other words, the FFLO 
state in the optical lattice does not have an internal de- 
gree of freedom in contrast to the continuum gases. It is 
interesting to study the effect of broken rotation symme- 
try by tuning the aspect ratio 3J]. Then, other FFLO 
states may be stabilized by the rotation. 

Third, the rotating FFLO phases studied in this pa- 
per are regarded as an analog of the superconducting 
states in the non-centrosymmetric systems [55i] . The de- 
generacy of the FF phases is lifted in the latter case by 
the anti-symmetric spin-orbit coupling arising from the 
broken inversion symmetry. The helical superconducting 
state is an analog of the FF state in our case, while the 
stripe phase corresponds to the intermediate state. How- 
ever, these non-centrosymmetric superconducting phases 
cannot be realized because the orbital de-pairing effect 
arising from the magnetic field affects the superconduct- 
ing state. In contrast, cold Fermi gases are ideal system 
for the study of FFLO phases thanks to their high con- 
trollability. 

A ring trap has been realized for Bose gases [H, HI] , 
and produces many intriguing topics. We expect that the 
realization of Fermi superfiuid on the ring trap will open 
a new research field of cold atom gases, as the devel- 
opment of sophisticated trap potentials triggered novel 
ideas. This paper shows that the rotating FFLO super- 
fiuid phases in the imbalanced Fermi gases will be one of 
the goals. We hope that active investigations will eluci- 
date such novel phases. 
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Appendix A: Derivation of the Hamiltonian 

In Appendix A, we derive the lattice Hamiltonian in 
Eq. ([ij from the continuum model. The Hamiltonian 
of continuum one-dimensional gases is described in the 
rotating coordinate as 



1 9^ 
'2M dx^ 

^\U\ j dx^{x)^l{x)i;^(x)^;^{x), 

Ho - \U\ f dx^{x)i^l{x)^^{x)^^{x), (Al) 



where i/^l-ix) and i/^aix) are the fermionic creation and 
annihilation operators, respectively, and Hq is the single- 
particle Hamiltonian. 

The coordinate along the ring x is given in term of 
the angle 9 and the radius of ring R. Then, the z- 
component of the angular momentum is described as 

With use of the Fourier transformation, the single- 
particle Hamiltonian Hq is obtained as 



k.CT 



(A2) 



Thus, the single particle energy is obtained as e{k) = 
— Rflk. It is understood that the Fermi momentum 
is shifted to the positive direction of fc for $7 > 0. On 
the other hand, the single particle energy in the lattice 
Hamiltonian (Eq. (H}) is obtained as 

e(/c) = -2tcos(fc - $). (A3) 

This is approximated in the continuum limit fc, $ <C 1 as 

e(fc) ~t($2 _2) _2t$A:-|-tA:2, (A4) 

Hence, we obtain the parameters of the lattice Hamilto- 
nian Eq. ^ by which the continuum gases are approxi- 
mately modeled. 



t ^ — , 2t$ o Rfl. 
2A/' 



(A5) 



Appendix B: Periodicity with respect to $ 



With use of the gauge transformation. 



the hopping term is transformed as 



— fr' r ■ 



'te^'^c' c- 



(Bl) 



(B2) 



Since we adopt the periodic boundary condition, the 
phase $ should be 



27r 



(B3) 



where n is an integer. This gauge transformation changes 
the Hamiltonian for the gas at rest to that for the rotating 
gas. This means that the physical quantities are periodic 
in $ with the period 2'k/Ni^ = $0- Hence we show the 
results for < <i>/$o < 1 in Figs. [2] and [T] 
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